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Understanding the physical origin of creep in calcium-silicate-hydrate (C-S-H) is of primary 
importance, both for fundamental and practical interest. Here, we present a new method, based 
on molecular dynamics simulation, allowing us to simulate the long-term visco-elastic deformations 
of C-S-H. Under a given shear stress, C-S-H features a gradually increasing shear strain, which 
follows a logarithmic law. The computed creep modulus is found to be independent of the shear 
stress applied and is in excellent agreement with nanoindentation measurements, as extrapolated to 
zero porosity. 


INTRODUCTION 

Creep is a major limitation of concrete. Indeed, it has 
been suggested that creep deformations are logarithmic, 
that is, virtually infinite and without asymptotic bound, 
which raises safety issues [1]. The creep of concrete is 
generally thought to be mainly caused by the viscoelas¬ 
tic and viscoplastic behavior of the cement hydrates [5]. 
While secondary cementitious phases can show viscoelas¬ 
tic behavior j2| , the rate and extent of viscoelastic defor¬ 
mations of such phases is far less significant than that 
calcium-silicate-hydrate (C-S-H), the binding phase of 
the cement paste |5]. As such, understanding the physical 
mechanism of the creep of C-S-H is of primary impor¬ 
tance. 

Despite the prevalence of concrete in the built environ¬ 
ment, the molecular structure of C-S-H has just recently 
been proposed HIS], which makes it possible to investi¬ 
gate its mechanical properties at the atomic scale. Here, 
relying on the newly available model, we present a new 
methodology allowing us to simulate the long-term creep 
deformation of bulk C-S-H (at zero porosity, i.e., at the 
scale of the grains). Results show an excellent agreement 
with nanoindentation measurements [5]. 

SIMULATION DETAILS 

To describe the disordered molecular structure of C-S- 
H, Pellenq et al. |3] proposed a realistic model for C-S-H 
with the stoichiometry of (Ca0)i.65(Si02)(H20)i,73. We 
generated the C-S-H model by introducing defects in an 
11 A tobermorite j7j configuration, following a combina¬ 
torial procedure. Whereas the Ca/Si ratio in 11 A tober¬ 


morite is 1, this ratio is increased to 1.71 in the present 
C-S-H model, through randomly introducing defects in 
the silicate chains, which provides sites for adsorption of 
extra water molecules. The ReaxFF potential [8], a reac¬ 
tive potential, was then used to account for the reaction 
of the interlayer water with the defective calcium-silicate 
sheets. More details on the preparation of the model and 
its experimental validation can be found in Ref. [S] and 
in previous works jH [SHU] ■ 

We simulated the previously presented C-S-H model, 
made of 501 atoms, by molecular dynamics (MD) using 
the LAMMPS package |17| . To this end, we used the 
REAXFF potential [Sj with a time step of 0.25fs. Prior 
to the application of any stress, the system is fully relaxed 
to zero pressure at 300K. 

RESULTS 

The relaxation of C-S-H, or of other silicate materials, 
takes place over long periods of time (years), which pre¬ 
vents the use of traditional MD simulations, which are 
limited to a few nanoseconds. To study the long-term 
deformations of C-S-H, we applied a method that has 
recently been introduced to study the relaxation of sili¬ 
cate glasses HH). In this method, starting from an initial 
atomic configuration of glass, formed by rapid cooling 
from the liquid state, the system is subjected to small, 
cyclic perturbations of shear stress ±Tp around zero pres¬ 
sure. For each stress, a minimization of the energy is 
performed, with the system having the ability to deform 
(shape and volume) in order to reach the target stress. 
These small perturbations of stress deform the energy 
landscape of the glass, allowing the system to jump over 
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FIG. 1. Shear strain and potential energy with respect to 
the number of loading/unloading cycles. The inset shows the 
shape of the applied shear stress. 


energy barriers. Note that the observed relaxation does 
not depend on the choice of Tp, provided that this stress 
remains sub-yield [TB]. This method mimics the artificial 
aging observed in granular materials subjected to vibra¬ 
tions |19| . 

Here, in order to study creep deformation, we add to 
the previous method a constant shear stress tq, such that 
Tp < To (see the inset of Figure [^. When subjected to 
shear stresses of different intensities, C-S-H presents a 
shear strain 7 that: ( 1 ) increases logarithmically with the 
number of cycles N (Figure]^ and (2) is proportional to 
the applied shear stress (see Figure]^. 

DISCUSSION 

The creep of bulk C-S-H can then be described by 
a simple logarithmic law 7 (iV) = (To/C)log(l -|- N/Nq), 
where Nq is a constant analogous to a relaxation time 
and C is the creep modulus. A careful look at the inter¬ 
nal energy shows that the height of the energy barriers, 
through which the system transits across each cycle, re¬ 
mains roughly constant over successive cycles. Accord¬ 
ing to transition state theory, which states that the time 
needed for a system to jump over an energy barrier Ej\ is 
proportional to exp(—we can assume that each 
cycle corresponds to a constant duration At, so that a 
fictitious time can be defined as t = NAt m- 

We note that the computed creep moduli C does not 
show any significant change with respect to the applied 
stress tq. As such, it appears to be a material prop¬ 
erty that can directly been compared to nanoindentation 
results extrapolated to zero porosity [ 6 ]. As shown in 
the inset of Figure we observe an excellent agreement, 
which suggests that the present method offers a realistic 
description of the creep of C-S-H at the atomic scale. 



FIG. 2. Shear strain with respect to the number of load¬ 
ing/unloading cycles, under a constant shear stress of 1, 2, 
and 3 GPa. The inset shows the creep modulus C with re¬ 
spect to the packing fraction (j> obtained from nanoindentation 
[0], compared with the computed value at (p = 0. 


This result also suggests that, within the linear regime 
(i.e., for sub-yield stresses, when C remain constant), de¬ 
formations due to cyclic creep and basic creep, with re¬ 
spect to the number of stress cycle or the elapsed time, 
respectively, should be equivalent. 


CONCLUSION 

We reported a new methodology based on atomistic 
simulation, allowing us to successfully observe long-term 
creep deformations of C-S-H. Creep deformations are 
found to be logarithmic and proportional to the applied 
shear stress. The computed creep modulus shows an ex¬ 
cellent agreement with nanoindentation data, which sug¬ 
gests that the present methodology could be used as a 
predictive tool to study the creep deformations of alter¬ 
native binders. 
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